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We present a mathematical model of glucose-lactose diauxic growth in Escherichia coli including 
both the positive and negative regulation mechanisms of the lactose operon as well as the inducer 
exclusion. To validate this model, we first calculated the time evolution of /3galactosidase for only the 
lactose nutrient and compared the numerical results with experimental data. Second, we compared 
the calculated cell biomass of the glucose-lactose diauxic growth with the experimental optical 
density of the diauxic growth for a particular E. coli MG 1655. For both cases, the numerical 
calculations from this model are in good agreement with these two experiments' data. The diauxic 
growth pattern of a wild type E. coli was also investigated. 



I. INTRODUCTION 

The complexity of how cells respond to their chang- 
ing environments requires the understanding of gene reg- 
ulation and protein expression. An interesting diauxic 
growth pattern of Escherichia coli (E. coli) on glucose 
and lactose nutrients has been extensively studied and 
well characterized for a long time (e.g., Jacob and Monod 

19 1, Jacob and Monod [2Cj], Beckwith [l[, Nester et al. 

23||). Fig. Q] illustrates a biphasic growth pattern of 
E. coli in a glucose-rich and lactose-rich medium. Since 
extra energy is required to make additional enzymes (lac- 
tose permease and /3galactosidase) for the lactose trans- 
portation and metabolization, the utilization of glucose 
was always preferred. The diauxic growth pattern can be 
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FIG. 1: Diauxic Growth of E. coli on glucose and lactose. I: 
glucose growth phase; II: lag phase; III: lactose growth phase. 
E. coli preferentially takes glucose and when the glucose nu- 
trient is depleted, a short lag phase appears before the cell 
growth resuming on the lactose utilization. Note that the 
logarithmic scale used in the number of cells. 
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explained by the lactose operon genetic regulatory net- 
work and the inducer exclusion. 

Quantitative mathematical models of the lactose 
operon and diauxic growth have been investigated by 
many authors (Jong [8(, Yagil and Yagil [32|, Wong et 
al. [3l|, Kremling and Gilles [llj, Kremling et al. 
Mahaffy and Savev [18| , Yildirim and Mackey 13311 . Satil- 
lan and Mackey (28[, Bailey 0, Lee and Bailev|lq| . Lee 
and Bailey [l6[ , Laffend and Shuler [l3[ , Ozbudak et al. 

A recent genetically structured model for the di- 
auxic growth of E. coli on glucose and lactose was de- 
veloped by Wong et al. [3lJ. In their work, all known 
information including the lactose operon transcription 
and translation processes, the flux models of unidirec- 
tional glucose transport and bidirectional lactose trans- 
port, the inducer exclusion, the catabolite repression, and 
the repressor DNA loop repression, was described by a set 
of nonlinear ordinary differential equations (ODEs). By 
solving the system ODEs with appropriate experimen- 
tal parameters, they produced a typical diauxic growth 
curve. However, no comparison was made with exper- 
imental data. Kremling and Gilles [ll[ and Kremling 
et al. pjj developed another structured mathematical 
model of glucose-lactose diauxic growth and it was vali- 
dated by their own experiments. 

Many other studies of lactose operon modeling were 
on the stability analysis of steady states and dynamic 
evolutions of E. coli on the lactose nutrient only. Ma- 
haffy and Savev [l8| and Yildirim and Mackey [33| de- 
veloped their mathematical models including the time 
delay of translation processes explicitly. Similar to Yagil 
and Yagil [13] , a simple Hill's function was applied as the 
regulation function of the lactose operon. Mahaffy and 
Savev [l8[ reported that there was a hysteresis of solution 
within a certain range of the parameter space. Yildirim 
and Mackey [33| showed a bistable steady state of the lac- 
tose operon. However, a common problem in their model 
was the assumption of a constant growth rate. This as- 
sumption is only valid in the exponential growth period, 
and may not be true through the period from the initial 
growth of E. coli to the exponential balanced growth. 

In this paper, we extended the mathematical model 
proposed by Wong et al. [3l| on one hand by includ- 
ing more realistic complexities (for example, the time 



2 



delay of translation processes of protein synthesis); on 
the other hand, we simplified their model (for example, 
by using an intuitive definition of the growth rate (Krem- 
ling et al. [12]). As is known, validation and prediction 
are two essential features for any useful mathematical 
model for real biological systems. Therefore, the focus 
of this present work is to validate this extended mathe- 
matical model by two particular experiments' data, and 
predict the diauxic growth pattern of a wild type E. coli 
on glucose and lactose in general. 



To construct the mathematical model, we made four 
changes to the model of Wong et al. [__tj]. First, a modified 
bidirectional flux model based on a simple microscopic 
description of the symmetric lactose transport was de- 
rived. This change connects the difference of the influx 
and efflux of the lactose transport. Second, the time de- 
lay of translation processes was considered explicitly. As 
pointed out by Yildirim and Mackey ,33j, it was impor- 
tant but failed to appear in Wong's model. In addition, 
the rate of translation was modeled by an empirical func- 
tion of the E. coli growth rate (Bremer and Dennis [_J]). 
The third change concerned the growth rate of E. coli. In 
our treatment, we used the growth rate that depends on 
the glucose and lactose uptake rate. This intuitive defi- 
nition let us eliminate the stiffness of the model of Wong 
et al. j[31J. The last change was on the inducer exclu- 
sion. A Hill's function with a higher steepness to model 
the inhibition of lactose permease was employed since it 
was reported by Kimata et al. [9( that the inducer exclu- 
sion may play a major role in the glucose-lactose diauxic 
growth pattern. 



RNAp RNAP 

Ijjo.) : . i : y 




Repressor Allolactose Pealactosidase Permease trausaeerylase 



FIG. 2: The lactose operon (Nester et al. [23,], Beckwith 0], 
Echols is a segment of DNA molecule in an E. coli cell 
with its length ~5000 base pairs (bps). It consists of three 
genes and three different regulation sites. The three genes are 
lacZ, lacY and lacA and their products are /3galactosidase, 
lactose permease, and transacetylase respectively. The en- 
zyme /3galactosidase breaks lactose and allolactose into glu- 
cose and galactose, and transforms lactose to allolactose. Lac- 
tose permease is in charge of transporting lactose into or out 
of the cell. The function of transacetylase is not clear in the 
lactose operon regulatory network, and thus is not considered 
in this study. Three different regulation sites are cyclic AMP 
(cAMP) receptor protein (CRP) binding site, promoter, and 
operator. These three binding sites provide physical locations 
for three transcription factors: CRP, RNA polymerase, and 
repressor protein. Upstream of the lactose operon, a repressor 
gene lad constitutively expresses the repressor protein. 



TABLE I: Effectors, transcription factors and binding sites in 
the lactose operon system 

Effector Transcription factor Regulation site 

a RNA polymerase promoter 

cAMP CRP CRP binding site 

allolactose repressor operator 



The studies presented here had three parts. In the 
first part, we defined three levels of /Jgalactosidase: 
basal level, partly induced level, and fully induced level. 
For the lactose nutrient only, we calculated the time 
evolution of /3galactosidase and compared it to Knorre 
[Io| experimental data. In the second part, we studied 
the diauxic growth of a particular E. coli MG 1655. 
The cell biomass of E. coli MG 1655 was calculated by 
inputting the same concentrations of glucose and lactose 
as those in experimental studies by Chang et al. [_| 
(http:/ /www. ou.edu/microarray/oumcf/growth_trans. htm). 
To fit the experimental data, we converted the cell 
biomass to the corresponding experimental optical den- 
sity (OD), and provided an analysis of the discrepancy 
between the numerical results and the observed OD 
in the lactose growth period. The numerical results 
of the mathematical model for both cases are in good 
agreement with these two experiments' data. In the 
third part, by adjusting the rate of glucose transport, 
the diauxic growth pattern of a wild type E. coli was 
studied. 



II. THE MODEL 

A. Biological model 

The structure of the lactose operon system is shown in 
Fig. [2] The transcription of the lactose operon is reg- 
ulated by the interaction of three transcription factors 
with their corresponding effectors at their binding sites 
in DNA listed in Table __ The a effector, if associated 
with the RNA polymerase, can help the RNA polymerase 
to find its promoter site and prepare to initiate the tran- 
scription. 

The lactose operon genetic regulatory network has 
both positive and negative feedback loops. The catabo- 
lite repression is positively regulated by the extracellular 
glucose. The high level of the extracellular glucose re- 
sults in a low level of cyclic adenosine monophosphate 
(cAMP), so cAMP receptor protein (CRP) is not acti- 
vated. If the extracellular glucose density is very low, 
then a high level of cAMP will result in an activation 
of CRP. cAMP as a modulator binds CRP to positively 
control its DNA binding activation. The complex of 
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FIG. 3: The metabolic pathway of glucose-lactose diauxic 
growth in E. coh 



tose, a plateau with a negligible growth rate shows up 
before the cell growth on the lactose utilization. The 
plateau shown in Fig. Q] characterizes the glucose-lactose 
diauxic growth in E. coli. 



B. Mathematical model 

The dynamics of the glucose-lactose diauxic growth 
has been modelled by Wong et al. jlH by using a set 
of nonlinear ODEs. Here, we present the four changes 
in the extended mathematical model which are differ- 
ent from the Wong et al. [3l| study. The complete rate 
equations to solve the dynamics of the glucose-lactose 
diauxic growth are listed in Appendix A. The units of 
the intracellular quantities and the extracellular quanti- 
ties are mol/gDC'W (gDCW: gram Dry Cell Weight) 
and M respectively. The units of the intracellular quan- 
tities can be converted to M by using the cell density 
Pceii — 300 gDCW I L when necessary. The unit of 
biomass is gDCW . 



cAMP~CRP at its binding site facilitates RNA poly- 
merase binding at the promoter site and helps RNA poly- 
merase in the transcription of the lactose mRNA. 

The negative regulation feedback loop is controlled 
by the interaction between the complex of repressor- 
operator and allolactose. If the intracellular allolactose 
(the nature inducer) is not available, the repressor pro- 
tein is bound at the operator site tightly, and thus the 
transcription of RNA polymerase is blocked. If the level 
of the intracellular allolactose is high, the allolactose will 
bind with the repressor and free the operator so that the 
transcription of the lactose o per on can be initiated. It 
was reported by Oehler et al. [25| that the binding affin- 
ity of the complex of repressor protein and allolactose at 
the operator site was lowered by about a few thousand 
fold. If the transcription of the lactose operon is dere- 
pressed, then /3galactosidase and lactose permease will 
be expressed. 

The glucose-lactose diauxic growth can be explained 
by the positive and negative regulation feedback loops 
and the inducer exclusion. The mechanism of the lac- 
tose transporter protein (lactose permease) inhibited by 
the interaction between lactose permease and the un- 
phosphorylated enzyme type II A subunit (EIIA), a 
component of phosphoenolpyruvate:sugar phosphotrans- 
ferase system (PTS) (Postma et al. [27]), is called the 
inducer exclusion, or glucose effects. Fig. [3] illustrates 
the metabolic pathway of E. coli in a glucose and lac- 
tose medium. When the extracellular glucose is negligi- 
ble, the inducer exclusion and the positive regulation are 
minimized. Because it takes time to induce the lactose 
operon to make the necessary enzymes /3galactosidase 
and lactose permease to transport and metabolize lac- 



1. Flux model of lactose transport 

The flux of the lactose transport entering the cell is 
given by, 



dt 



— k lac (p ce iip perm (t))—^- 



Plactt) ~ PcellpLc 



Placid + PceliPi (*) + Plac 
(1) 

where p° ac , p l lacl and p per m denote the densities of the ex- 
tracellular lactose, the intracellular lactose, and the intra- 
cellular lactose permease respectively. The rate constant 
of the lactose transport is kj ac and the saturation con- 
stant of lactose transport is p* ac . The flux of the lactose 
transport is proportional to the the density of lactose per- 
mease regulated by the lactose operon. The derivation 
of Eq. [T] based on the microscopic description is given in 
Appendix B. 

Eq. [1] is different from the flux model of the lactose 
transport given by Wong et al. [3l|. First, this flux model 
is proportional to the difference of the extracellular and 
intracellular lactose densities. Second, there is only one 
saturation constant for both influx and efflux of the lac- 
tose transport. In our numerical results shown later, the 
intracellular lactose density is very small compared to the 
extracellular lactose density. This could explain why the 
saturation constant of the efflux is very high in the efflux 
model of the lactose transport by Wong et al. [3l| . 



2. Time delay of translation 

Proteins such as /3galactosidase and lactose permease 
are synthesized in the translation process. The rate equa- 
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tions of them arc, 



dt 

dpperm (^) 

dt 



-k t sPrnRNA{t ~ Tp) - (p + d p ) p p ga l{t\2) 



ktsPmRNA(t 



(p + d p )p perm (t)(3) 



where p m RNA denotes the density of mRNA molecules. 
rp and t p are the time delay constants of /3galactosidase 
and lactose permease respectively. These two constants 
are the results due to the finite speed of producing 
polypeptides at the ribosome sites. The factor \ in Eq. [5] 
accounts for a /3galactosidase molecule consisting of four 
monomers. The number of proteins produced per unit 
time, i.e., the rate of translation constant, is kt s . This 
rate is dependent on the growth rate ([H, [IB]). Here, the 
relationship between k ts and the growth rate pit) was 
assumed to be, 
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where is the growth rate corresponding to a one-hour 
doubling time. If there is no lactose in the nutrient pool, 
then the lactose operon is maintained at its basal level by 
the stochastic and rare transcription. Thus, the rate of 
translation k ts was assumed to be k^ sal . If lactose nutri- 
ent is available, then the rate of translation is increased 
until p > p . The higher the steepness parameter m is, 
the faster the rate of translation k ts returns to its nor- 
mal level from k\^ sal . The protein degradation time of 
/3galactosidase and lactose permease is d p . 



3. Inducer Exclusion 

When a glucose molecule is transported into a cell, it 
is phosphorylated into a glucose-6-phosphate by a phos- 
phorylated EIIA. The unphosphorylated EIIA inhibits 
the lactose permease and thus prevents the lactose trans- 
port. This inhibition is currently known only for the ex- 
tracellular lactose. Therefore, in the expression of the 
flux model of the lactose transport, an inhibition factor 
is attached to the influx of the extracellular lactose den- 
sity. We used a similar Hill's function (Wong et al. [3l| ) 
but with a higher steepness for modeling the inducer ex- 
clusion process in this study. It was reported by Kimata 
et al. that glucose effects may play an important role 
in the diauxic growth pattern. To incorporate this infor- 
mation, the inducer expression is modeled by, 
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where the inhibition constant is Ki ex and the extracel- 
lular glucose density is p° g i u - The total flux of lactose 
coming into the cell with the inducer exclusion is now 



modified, 
dt 



^lac [PcellPperm 

(*))^ 



PlacitfRiex ~ PcellPlac 



Plac^f) + PcellP] (*) + 

Plac 

(4) 



4- Growth Rate 



The doubling time may be defined as the time needed 
for a cell absorbing the number of glucose molecules 
( N glu) t0 make a celL To nnd this number N c g f u , the 
yield coefficient of glucose (Yx/giu) an d an average E. coli 
cell biomass were used. The doubling time is thus given 
as a function of the flux of glucose or lactose transport, 



T,, 



N cell / dNgl u 

9lu/ dt • 



and the growth rate is defined as p = fe^. The intra- 
cellular lactose and allolactose is hydrolyzed into glucose 
and galactose, and galactose is then converted into glu- 
cose. We simply used one lactose molecule contributing 
two glucose molecules. The growth rate is thus given, 



p(t) = 



fag ,. t PghM 



glu 



PglJt) + PW 



n U\ RiexP° ac (t) P\ ac {t) . , , 

+ 2k lacPPgal{t) - j Ts ■ 77s —)■ (5) 

PlacW + PlacW + Plac 

The rate of the glucose transport and the saturation con- 
stant of the glucose transport are denoted as k glu and 
Pgiu respectively. Note that the unit of N' g f^ 11 is now 
mol/gDCW. The expression of the growth rate shows 
that p{t) depends on the extracellular glucose and lactose 
densities as well as the inducer exclusion. This intuitive 
definition of the growth rate (Kremling et al. [l2j]) avoids 
the stiffness resulting from the assumption of the growth 
rate depending on the level of glucose-6-phosphate and 
the extra burden of the lactose operon (Wong et al. 



III. RESULTS 

To solve the dynamics of E. coli growth on glucose and 
lactose nutrients, the rate equations in Appendix A were 
integrated by the fourth-order Runge-Kutta algorithm. 
All the parameters used in the rate equations are listed 
in Table HU 

We first assumed the volume of the glucose or lactose 
nutrient pool to be infinite but their densities were kept 
as constants all the time during the cell growth. This 
corresponds to the average cell growth studies without 
considering the change of the extracellular glucose or lac- 
tose densities. Second, we took the volume of glucose or 
lactose pool as 1 I, but the extracellular glucose and lac- 
tose densities were varied with time. The dynamics of the 
glucose-lactose diauxic growth were then investigated. In 
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FIG. 4: Calculated curves (solid line) from the ex- 
tended mathematical model, from top to bottom, p° ac = 
40, 4, 2, ImM, were compared to Knorre's experimental data 
(open circle). 



both numerical calculations, the repressor density was 
assumed to be a constant, which is roughly 15 repressor 
proteins per cell. Finally, we studied the diauxic growth 
pattern of a wild type E. coli by adjusting the rate con- 
stant of the glucose transport. 



A. Lactose growth 

The induction status of the lactose operon in this study 
is described by the basal level, partly induced level, and 
fully or completely induced level. The fully induced sta- 
tus is defined as the maximal level of /3galactosidase when 
the lactose operon works at its maximal transcription ef- 
ficiency. It was reported that the repression level of the 
lactose operon increased 2400 - 3000 fold for 10 - 15 re- 
pressor proteins per cell ([13]). Therefore, we defined 
the basal level of pj^f as about one 2700th of pfig%. 
The maximal value of /3galactosidase was obtained by 
letting the extracellular lactose density be very high. A 
partly induced lactose operon was defined as being be- 
tween the observable induction level and the maximal 
induction level of /3galactosidase. Similar to Ozbudak 
et al. [Ill, this observable induction level was defined as 
one 100th of p™g£v ^ the induction level of the lactose 
operon falls between the basal level and the observable 
level, the lactose operon is not induced. If the induction 
level is above the observable level but less than the max- 
imal levels, then the lactose operon is partly induced. 

In the numerical calculations, the maximal level of 
/3galactosidase, lactose permease, intracellular lactose 
and allolactose were obtained by letting p° ac = 4 x 
10~ 2 M or above, because the maximal transcription 
efficiency was saturated at the very high extracellular 
lactose density. By using the parameters listed in Ta- 



ble OH pf<% is 2.1 x 10~ 8 mol/gDCW, which is about 
3150 /3gafactosidase molecules per cell if 1 molecule per 
cell corresponds to 2 x 10~ 9 M. First, we let the ex- 
tracellular lactose density be zero, and solved the sys- 
tem ODEs. The basal level of the rate of translation 
fcbasai wag ac y us t c d anc [ estimated such that the induc- 
tion level of the lactose operon was less than or equal to 
the basal level. Next, a few numerical experiments were 
done within a range of the extracellular lactose densities 
from 1 x 10 -3 M to 4 x 10~ 2 M. The calculated results 
of /3galactosidase were normalized by their steady state 
values respectively. In Fig. [H the numerical curves of 
/3galactosidase for four different levels of the extracellu- 
lar lactose were plotted and compared to Knorre's exper- 
imental data (PIl)- At the very beginning (less than 25 
minutes) of the cell growth on lactose only, the sigmoid 
shape of these calculated curves can be seen in Fig. [4] 
This is in contrast to the hyperbolic sha pe g iven in the 
previous study by Satillan and Mackey [28(. The rea- 
son for this change is mainly due to the change we made 
on the growth rate /x, i.e., p(t) is the time-dependent 
growth rate while it was assumed to be a constant in 
the previous study (Satillan and Mackey (28J). Other 
changes such as the time delay of the protein transla- 
tion and the time-dependent translation rate (depending 
on the growth rate p(t)) also played some roles. The 
sigmoid shape of /3galactosidase could partly explain the 
lag phase of the glucose-lactose diauxic growth before the 
cell growth switched on lactose. 

If both glucose and lactose are present in an infinite 
pool, then the extracellular lactose and glucose densities 
are to be assumed to be fixed all the time during E. coli 
growth. In this situation with the initial cell biomass or 
cell number assumed to be very small, bistability of the 
lactose operon can be studied. A numerical study of the 
lactose operon bistability with full consideration of inter- 
actions of regulation sites on DNA with their effectors 
and transcription factors has been given by Satillan and 
Mackey [2^|. In their study, the level of /3galactosidase 
was taken as the order parameter to indicate the lac- 
tose operon transition from the uninduced state to the 
induced state. When the extracellular lactose density 
increased up to a threshold value, the level of pp ga i in- 
creased discontinuously from zero to a finite value. When 
the process was reversed, the level of /3galactosidase de- 
creased continuously through that threshold value. 

In this study, we did the numerical calculations on 
both increasing and decreasing extracellular lactose den- 
sity. However, we did not find discontinuation of Pf3 ga i- 
A possible reason of failing to see the bistability is due 
to the growth rate p. The growth rate in our study is 
time dependent and less than the constant growth rate 
used in Satillan and Mackey [H)]. Moreover, the tran- 
scription efficiency of the lactose operon in their calcu- 
lation was carried out by the partition function which 
was different from our approach. For lactose only in a 
nutrient pool, we found the threshold value to turn on 
the lactose operon at the observable induction level to be 
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FIG. 5: Comparison of the calculated optical density (solid 
line) with the experimental optical density of E. coli MG 1655 
(open circles). 



~ 3.6 fiM. If the extracellular lactose density is lower 
than this value, then the lactose operon is not induced; if 
the lactose density is higher than this value, then the lac- 
tose operon is partly induced. The fully induced lactose 
operon can only be reached by a high level of the extra- 
cellular lactose density. Ozbudak et al. [26[ studied the 
bistability of the lactose operon when both glucose and 
thio-methylgalactoside (TMG) were used as substrates. 
The threshold value to turn on the lactose operon was 
estimated as 3 fj,M with only TMG. 



B. Glucose- lactose diauxic growth 

If the glucose and lactose nutrient pool is finite, then 
the extracellular glucose and lactose densities will be 
changed during the growth period of E. coli. There have 
been many experimental works on the glucose-lactose di- 
auxic growth pattern, and most experimental data have 
been given in terms of OD rather than the cell biomass 
(e.g., [Hi Hi)- [H studied the glucose-lactose diauxic 
growth; the experimental OD was measured at 600 nm 
for the E. coli MG 1655. The MG 1655 was cultured 
in 1 liter with a medium containing 0.5 g/l glucose and 
1.5 g/l lactose. To monitor the cell growth, the OD at 
600 nm was measured. 

In numerical calculations, the glucose and lactose ex- 
tracellular densities (0.5 g/l and 1.5 g/l) correspond to 
Pglui* = 0) = 2.78 mM and p° la Jt = 0) = 4.38 mM, 
respectively. It was reported that the lag phase or the 
transition between the glucose growth and the lactose 
growth lasted about 36 minutes. The OD discussed by 
Lawrence and Maier 11411 is accurate when its values in- 



crease from 0.1 to 0.6, and is linearly related to the cell 
biomass. If OD is higher than 0.6, the true OD should 
be higher than that observed in the experimental mea- 
surement, and the linear relationship is no longer valid. 
The real OD can be corrected by a standard curve in 
a specific experiment. Since the rate of glucose uptake 
depends on the extracellular glucose density, the num- 
ber of EIIBC, and the dwelling time associated with the 
phosphorylation process, different species of E. coli may 
have different numbers of EIIBC and dwelling time. This 
results in the different growth rate and doubling time. 

The typical doubling time in glucose exponential 
growth of E. coli MG 1655 is about 60 minu tes 
( [http://www.genome.wisc.edu/resources/strains.html ). 
To fit the experimental data given by Chang et al. [H, 
the rate constant of glucose transport, k l glu , was taken as 
2.1 x 10~ 4 mol glu / gDCW / min while a typical range of 
this rate constant is 1.0 - 3.2x 10~ 4 mol glu / gDCW / min. 
The rate constant of glucose transport determines the 
doubling time of the cell growth on glucose (see Eq. 

EH). 

The calculated OD in Fig. [5] is in good agreement 
with the experimental OD during the glucose growth and 
the lag phase, but it is higher during the lactose growth. 
This discrepancy can be explained by the OD measure- 
ment mentioned before. Since the OD during the lac- 
tose growth is greater than 0.6, the real OD should be 
higher than the experimental measured OD (Lawrence 
and Maier (HI), and a standard calibration can be ap- 
plied to correct the measurement due to the insensitiv- 
ity of the light scattering. This statement actually con- 
firms that the higher calculated curve of OD is consistent 
with the analysis of the real OD measurement. Because 
the initial cell mass given in this study is very small, 
10~ 3 gDCW , we have shifted the time and used a pro- 
portional constant between the cell biomass and the OD. 
The proportional constant between the OD and the cell 
biomass is 2.62, i.e, OD = 2.62X. In Fig. [5]the lag phase 
was estimated to be 32 - 38 minutes. 



C. Diauxic growth pattern of a wild type E. coli 

The doubling time of the typical wild type E. coli is 
usually within 20 to 40 minutes, which is lower than the 
doubling time of E. coli MG 1655. To reflect this charac- 
teristic, we increased the rate constant of glucose trans- 
port (k glu = 3.2 x 10~ 4 mol glu / gDCW / min) for a wild 
type E. coli. 

In the numerical calculations, four levels of the initial 
extracellular glucose and lactose were set as: p° lu (t = 
°) = Piad* = 0) = 1, 2, 3, 4 mM. The initial levels 
of mRNA, /3galactosidase, and lactose permease were set 
to their basal levels. The biomass was initially used as 
10 -5 gDCW . All other quantities were zero initially. 

In Fig. \§\ we plotted time courses of cell biomass, 
growth rate, mRNA, /3galactosidase, intracellular lactose, 
extracellular lactose and glucose. The diauxic growth 
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patterns (logarithmic scale used in the cell biomass) on 
glucose and lactose are shown in Fig. [5] (a) . The slopes 
of the cell biomass represent different growth rates. 

From Fig. [S](b), during the exponential growth period 
of glucose, the growth rate was almost a constant. 
However, when the level of the extracellular glucose de- 
creased after 500 minutes, the growth rate /i(t) dropped 
rapidly. In the following lag phase, the growth rate was 
negligible or less than the threshold value of detecting cell 
growth. After the lag phase, E. coli restarted the lactose 
utilization and the growth rate increased rapidly. When 
the lactose nutrient was almost used up the growth rate 
dropped quickly again. 

In Fig. [6] (b), during the lactose growth period, the 
growth rate fi(t) never reached a constant, so no ex- 
ponential growth was reached. This is accounted for 
by the high cell biomass accumulated after the glucose 
growth period. The base cell biomass was now (5 x 10~ 2 
- 2 x 10 _1 gDCW. Therefore, during the lactose growth 
period, there were no enough lactose nutrient supplied 
before the lactose exponential growth was reached. 

Fig. [U (c) shows how protein /3galactosidase changed 
with time. The protein expression increased rapidly 
once the lactose utilization started. However, the 
/3galactosidase took a long time to decrease to zero even 
if there were no lactose nutrient available. This can be 
explained by the protein degradation time which is ap- 
proximately 100 minutes. The time evolution of lactose 
permease was the same as /3galactosidase except at the 
level four times greater than the level of /3galactosidase. 

The level of mRNA zya is shown in Fig. M (d). The 
mRNA 2l)a molecules increased and decreased similar to 
the growth rate /j(t) during the lactose growth period. 
The peak densities of mRNA zya in Fig. [6] (d) were about 
15 - 18 molecules per cell. The level of the intracellular 
lactose in Fig. [6] was very low and it had the same time 
evolution as mRNA Z j, a . The time courses of the extra- 
cellular glucose and lactose densities during the glucose- 
lactose diauxic growth are plotted in Fig. [S] (f). From 
this figure, the extracellular lactose was not transported 
into the cell before the glucose was used up, and this can 
be attributed to the inducer exclusion. 



IV. SUMMARY 

We have presented an extended mathematical model of 
the glucose-lactose diauxic growth in E. coli. The math- 
ematical model made with four changes was validated by 
two experiments' data. For the lactose nutrient only in 
an infinite pool, the numerical result of time evolution of 
/3galactosidase was compared with Knorre's experimen- 
tal data. For both glucose and lactose nutrients in an 
finite pool, the dynamics of the cell biomass were con- 
verted to the experimental OD of E. coli MG 1655. The 
discrepancy between the calculated OD and the observed 
experimental OD in the period of the lactose growth was 
analyzed. The diauxic growth pattern of a wild type of E. 



FIG. 6: Diauxic growth pattern of wild type E. coli with four 
levels of the extracellular glucose and lactose (p° g i u {t = 0) = 
P°a C {t = 0) = 1, 2, 3, 4 mM). 
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coli was studied. The extended mathematical model in 
this work shows that diauxic growth patterns of different 
species of E. coli can be quantitatively investigated and 
understood well by integrating all available information 
of the involved biological processes with different time 
scales. 



dp° glu (t) 
dt 

dpf ac (t) 
dt 

dt 

dpi3gal{t) 

dt 

dp perm \p) 

dt 

dp\ac(t) 

dt 



dpallojt) 

dt 

dpcAMP (t) 

dt 

dX(t) 
dt 



Eq. [AT] describes the flux of glucose transport entering 
the cell. The extracellular glucose density is p° lu (t) . The 
rate of glucose transport and the saturation constant are 
Kin and P*giu respectively. 

Eq. IA2I describes the flux of lactose transport entering 
the cell. The extracellular lactose density is p° ' (t), and 
the intracellular lactose density is p\ ac {t). The cell den- 
sity p ce u is used for converting the unit (mol/gDCW) 
of the intracellular quantities having the same unit (M) 
as the extracellular lactose density. The rate of lactose 
transport and the saturation constant are k\ ac and p\ ac 
respectively. 

Eq. IA3I describes the transcription process of the lac- 
tose operon. The densities of mRNA and the lactose 
genes are p m RNA(t) and piaczya- The degradation rate 
and the rate of transcription are d m RNA and k tc respec- 
tively. The maximal transcription rate k tc is estimated as 
2/s, for the size of mRNA polymerase with the a factor 



APPENDIX A: RATE EQUATIONS 

The nonlinear ODEs (Wong et al. [3l[) used to solve 
the dynamics of the diauxic growth with the four changes 
made in this study are: 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 



(A6) 
(A7) 



(A8) 
(A9) 
I 

is ~ 37 nm and the elongation rate of mRNA molecule 
is about 50 bps/s. 

The factor P, representing the transcription efficiency 
of the lactose operon, is decomposed into P = Pi x P 2 x 
P3. Pi is the probability of the complex of RNA poly- 
merase with the a factor at the promoter site; P2 is the 
probability of the complex cAMP and CRP at the CRP 
binding site; and P3 is the probability of the empty op- 
erator site. The calculations of Pi, P2 and P3 are given 
in Appendix C. 

Eqs. IA4I and IA5I describe the rate changes of 
/3galactosidase and lactose permease. The rate of trans- 
lation is k ts and it depends on the growth rate p(t). 
The degradation rate is d p . Time delay constants of 
/3galactosidase and lactose permease are 7/3 and t p re- 
spectively. 

Eqs. IA6I and IA7I describe the rate changes of the in- 
tracellular lactose and allolactose. The hydrolysis of lac- 



PcellKgl u — 



' P°gl U (t) + P* 9 lu' 

_yt 1 r.\\ P°ac(t)Riex ~ Pcelip\ ac 

K lac\PcellPperm\l>)) ,,\ . j /,\ . * j 

PlacW + PcellPlacW + Plac 

hcP Piaczya ~ (m(*) + d m RN A ) PmRN A (t) , 

-rktsPmRNA{t - Tp) - {p(t) + dp)p/3 ga l(t) 

hsPmRN A(t — T p ) - (p(t) + d p )p perm (t), 
P° lac {t) ~ PceUP\ ac (t) 



klacPperm{t) Q 



-k?ac°PP g al{t) 



PUt) + PceliPlacit) + Pt 

pL(*) 



P\aM + K la 



-kfS*P0gal(t) , ~ P P \ ac {t\ 

Plac\ l ) ^'ac 

KuoP^it) -P(t)Pallo(t), 

PaUo\t) + K a u 

, K C AMP 
KcAMP 



K cAM p + P° glu [t) 

(k ex + p{t))p c AMp{t), 



p(t)X(t). 
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tose and allolactose was assumed to follow the Michaelis- 
Menten mechanism. The saturation constants of the lac- 
tose and allolactose hydrolyzation are Ki ac and K a u . 
The rate constants of transforming lactose or allolactose 
to glucose and galactose are kf£* and k^* respectively. 

Eq. IA8I describes the catabolite repression regulated 
by cAMP. The cAMP synthesis constant and inhibition 
constant are k c AMP and K c amp respectively. Eq. IA9I 
describes the dynamics of the cell biomass. The values 
of all parameters in the rate equations are given in Table 

IE 

Table HJ inserted here. 



APPENDIX B: DERIVATION OF THE FLUX 
MODEL OF THE LACTOSE TRANSPORT 

The microscopic description of bidirectional lactose 
transport is similar to that of the glucose transport; how- 
ever, the lactose transport is more complicated. For sim- 
plicity, we first derived the flux model of glucose trans- 
port and then the flux model of bidirectional lactose 
transport. 

Glucose is transported by the glucose transporter en- 
zyme type II BC (EIIBC), a component of PTS. In this 
study we assumed two states of EIIBC: either free or 
occupied. If the EIIBC subunit is occupied by a glu- 
cose molecule, then it can not accept any more glucose 
molecules until it is free again. The total number of EI- 
IBC on the plasma membrane of E. coli is given at time 
t, 

Nl t {t) = jv;(t) + Nl{t), 

where numbers of free and occupied by glucose are Nf (t) 
and Ng(t) respectively. The total flux of nutrients into 
the cell depends on the number of free EIIBC and the 
extracellular glucose density, 

Jcell = NJK c p° glu (t). 

K c is the proportional constant characterized by the na- 
ture of biochemical reactions of glucose and EIIBC. At 
any instant time, N%(t) glucose transporters are occu- 
pied by glucose molecules and if the time for a glucose 
molecule to be completely transported into the cell is t s , 
we then expect the number of glucose molecules per unit 
time entering the cell to be, 



The total glucose flux becomes, 



P — 

J cell — 



_ * 

t. 



The time t s for transporting a glucose molecule into cell 
may be associated with the phosphorylation process, and 
the time scale of this transportation process is ~ ms. 
Because it is very small compared to the time scale of 
changing the density of EIIBC, which is ~ s, these fluxes 
must be equal. This gives us the number of free EIIBC 



as, 



Nf(t) 



N t e ot (t) 



J, 



If the characteristic density p* lu is defined as, p* lu — 
l/(t s K c ) and l/t s — k g is identified as the rate of 
the complex of EIIBC and glucose releasing the glucose 
molecule into the cell, then the total glucose flux coming 
into an average cell is, 



Jceii = k g N t e ot (t) 



P° g iu(t) 



PgluM+Pglu 



This is the flux model usually used in modeling glucose 
transport. The characteristic glucose density actually is 
the saturation constant of glucose transport through EI- 
IBC. The equation of describing the rate of the extracel- 
lular glucose density is then, 



dp°i u (t) 



dt 



-Pcellk g i u 



P° g lu(t) 



PgmW+Plm 



where k glu is the rate constant of the glucose transport 
including the contributions of all glucose transporters on 
the E. coli cell plasma membrane. 

The transport of lactose is similar to the glucose trans- 
port. However, the lactose transport is bidirectional and 
cotransported with a proton ion through lactose perme- 
ase. Here the kinetic constant in the lactose transport 
was assumed to include the effect caused by the electric 
potential across the cell plasma membrane. The lactose 
permease was assumed to have three states: free (f), oc- 
cupied by lactose from outside the cell (oo) , and occupied 
by lactose from inside the cell (oi). Symmetrical lactose 
transport was also assumed. The total number of lactose 
permease in average cell is given as, 

NUt) = N l f (t) + N l 00 (t) + N l oi (t)- 

By using the same argument as in the glucose transport, 
the flux of lactose entering an average cell is now, 



Jcell = hN l tot {t) — 



PtacW ~ PcellPlac 



l + t s K c *p° lu (t) 



PlacW + PcellPlcW + Pic 

where p* ac is the characteristic constant or saturation 
constant of the lactose transport. 

The argument of two different time scales discussed 
above is in essence the equilibrium assumption of the glu- 
cose and lactose transportation at the longer time scale 
of changing the density of nutrient transporters. Indeed, 
the total number of glucose transporters is maintained at 
a constant number because the degradation is balanced 
with the expression of the glucose transporter. The num- 
ber of lactose transporters is regulated by the lactose 
operon so that the net number of lactose permeases is 
gradually increasing. Although the total number of the 
lactose permeases is not a constant at the longer time 
scale of changing the density of lactose permease, it can 
be viewed as a constant at the shorter time scale of trans- 
porting lactose into the cell. 
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APPENDIX C: CALCULATION OF 
TRANSCRIPTION EFFICIENCY 

The transcription initiated by mRNA polymerase will 
not start until the lactose operon is activated. The acti- 
vation of the lactose operon is mainly determined by the 
transcription efficiency (Bailey [2|). Presently, three ap- 
proaches are available. The first is to use a Hill's function 
as the regulation function and ignore the detailed molec- 
ular biochemical information of the lactose operon. This 
approach has been developed for a long time to simulate 
the regulation process of the lactose operon. The second 
approach is based on the assumption of fast equilibrium 
approximation; the transcription efficiency is treated as 
the product of three binding probabilities: the probabil- 
ity (Pi) of mRNA polymerase associated with a factor 
binding at the promoter site; the probability (P 2 ) of a 
complex of CRP and cAMP binding at the CRP site; 
and the probability (P3) of the free operator. From the 
mass balance law and given detailed biochemical reac- 
tions of the lactose operon, a set of nonlinear algebra 
equations can be constructed with all kinetic parameters 
from the experimental studies. The third approach is 
through identifying all possible states occurring in the 
binding sites and then giving a partition function con- 
sidering the statistical mechanics with the rapid equilib- 
rium approximation. The transcription efficiency, or the 
induction probability, is calculated by the states allowed 
for transcription to all possible states. 

The last two approaches are better than the first one 
because the genetic structured information of the lac- 
tose operon is incorporated. Therefore an accurate de- 
scription of the growth pattern based on the more de- 
tailed genotype information is expected. Nonetheless, a 
large set of kinetic parameters limits the application of 
the last two approaches for calculating the transcription 
efficiency. In this study, we used the second approach 
but only the main operator was considered. In addition, 
we assumed no dependence among these three processes. 
The total induction probability is thus, 

P = P 1 x P 2 x P 3 , 

where 

p(mRNAp,a,p) 
1 P(P) 

_ p{CRP,cAMP, CRP site) 
2 ~~ p{CRPsite) ' 

= p(o, free) 
pip, tot) 

Similar to Lee and Bailey [l5[ , Laffend and Shuler [l3[ , 
and Wong et al. [3l| , the transcription efficiency was cal- 
culated by three fractions. At any instant time, rapid 



equilibrium of the lactose operon was assumed. The bio- 
chemical reactions associated with this rapid equilibrium 
are written as, 

Ci 

Pxl + Px2 ■■— Pxlx2 

C 2 

Pxl + Px3 ^ Pxlx3 

C 3 

Pxl + Pxi ^ PxlA 

pxlx2 + Px3 ^ Pxlx2x3 
C5 

Pxlx2 + Pxi •<— Pxlx2x4- 

The equilibrium constants C\, C 2 , C3, C4, and C5 are 
listed in Table HTH The x\, X2, and x% are transcription 
factor, effector, and response element respectively (see 
Table [TJ) . X4 stands for the nonspecific DNA sites in E. 
coli. If the mass balance law is implied at any time, we 
have the following equations, 

Pxl [tot) = Pxl + Pxlx2 + Pxlx3 + Pxlx4 + Pxlx2x3 + Pxlx2xi 
Px2(tot) — 

Px2 ~t~ Pxlx2 Px\x2x3 ~t~ Pxlx2x4 
Px3{tot) = Px3 + Pxlx3 + Px\x2x3 
Pxiitot) — 

PxA "f" PxlxA "f" Pxlx2xA- 

Since the number of nonspecific sites of DNA is much 
greater than that of regulation sites, p X 4 is approximated 
as the total number of DNA sites, Pxiitot) = p"^™. The 
combinations of x\, Xi, X3, X4 represent the complexes 
of them. In this study, only complexes like x%X2, X1X3, 
X1X4, X1X2X3, and X1X2X4 exist. 
Table Hn] inserted here. 

By using equilibrium relationships, p x ix2 PxlPx2 j 

Pxlx3 = C2PxlPx3i Pxlxi = C^PxiPxi, p x \x2x3 = 

Ciip x ix2Px3, and = C^p x ix2Pxi, three nonlinear 

algebraic equations with three unknown quantities, p x i, 
Px2, and p x3 are given, 

Pxl(tot) = pxl + CiPxlPx2 + C 2 PxlPx3 + C 3 p x lPx4 
+C4C1 lPx2Px3 + C^C\PxlPx2Px1 

Pxiitot) = Px2 + CiPxlPx2 + C±CiPxlPx2Px3 
+ C§Cip x lPx2PxA 

Px3(tot) = Px3 + C 2 PxlPx3 + CiCiPxlPx2Px3 

The p x \{tot), p X 2{tot), and p X 3(tot) are input densities, 
which can be constant values (see table HU, or predicted 
by the differential equations [13j, whereas p x \, p X 2, and 
p X 3 are free parts of them. The transcription efficiency 
P = Pi x P 2 x P3 can be calculated by solving three sets of 
the above three nonlinear algebra equations respectively. 
For Pi, p Xl = p™ ap , p X2 = p a tot , and p Xa = (% ot \ for P 2 , 

crp /,\ j crpsite r j-> 

Pxi = Ptof Pxi = PcAMP(t), and p X3 = p to t ; tor P 3 , 

Pxi = Prepressor , Px 2 = Pallo(t)> and p X3 = Ptof TllC 

densities of cAMP and the intracellular allolactose are 
provided by the nonlinear ordinary differential equations 
in Appendix A. 
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TABLE II 


: Parameters used in modeling E. coli \ 


growth. 


parameters 


Values 


Itelerences 


it 

K glu 


1.0 - o.z X 1U mot glu/(gJJCW min) 


Tool 

[22| 


Pglu 


lo jllVl 


[27| 




1 A O ' „ — 1 

z!4o mm 


[31| 


Plac 


z.O X 1U JkZ 


[17| 


hi,; 


i on /mi'-n — l 
IzO vn%Ti 


estimated 


d m RNA 


u.oyo mm 


HO] 


L.eff 

"-ts 


lo.O TTim 


[lbj 


T P 


0.05 min 


Tool 
133] 


T p 


1.05 min 


Tool 
133] 


,] 

a p 


U.01 min 


15] 


7 basal 
K ts 


1.8 min 1 


estimated 


&lac 


1.4 X 1U IVl 


131] 


i^allo 


z.o X 1U iV'i 


ro -i i 
[3jJ 


h,cat 
^lac 


yo4U mm 


[31J 


■i alio 
K lac 


o4DU mm 


131] 


h.cat 
^allo 


i.o x iu min 






o v in - 5 ji/f 
z X 10 M 


[12| 


KcAMP 


1 w in — 5 71//" tan~ta<n 

i x iu iva /mm 


I'J 


KcAMP 


A v i n— 4 a 


\ r )A'\ 
124] 


h 

™- ' X 


2.1 min 1 


HI 


/ glu 


90 gDCW/(mol 


estimated 


Pcell 


300 gDCW/L 


[31] 


Placzya 


8.47 x 10" 12 mol/gDCW 


[31] 


Ptot 


8.47 x 10" 10 mol/gDCW 


estimated 


map 
Ptot 


5 x HT 6 M 


[31] 


Ptot 


1 x 10~ 6 M 


[31] 


crp 
Ptot 


2 x 10~ 6 M 


[31] 


crpsite 
Ptot 


8.47 x KT 12 mol/gDCW 


[31] 


o 

Ptot 


8.47 x 10" 12 mol/gDCW 


[31] 


non 
Ptot 


0.0118 M 


[31] 



TABLE III: Equilibrium constants involved in the biochemi- 
cal reactions of lac operon regulation sites ([3lJ, [HI). 



Parameters RNAP(:ei) 


^jixr{Xi ) 


Repressor ( x \ ) 




cr(x 2 ) 


cAMP(a; 2 ) 


Allolactose (x 2 ) 




Promoter^) 


CRP site (x 3 ) 


Operator(a;3) 


Ci 


2 x 10 9 M _1 


1 x 10 3 M _1 


2 x 10 6 M _1 


c 2 


1 x 10 6 M" 1 


1.1 x 10 5 M' 1 


1 x 10 13 M' 1 


c 3 


1 x 10 6 M' 1 


5 x 10 3 M" 1 


1 x 10 4 M" 1 




1 x 10 9 M' 1 


1 x 10 10 M" 1 


1 x 10 10 M- 1 


C 5 


1 x 10 5 M' 1 


1 x 10 5 M' 1 


1 x 10 s M _1 



